Microscopic Route to Nematicity in Sr 3 Ru 2 07 
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An anisotropic metallic phase dubbed electronic nematic phase bounded by two consecutive 
metamagnetic transitions has been reported in the bilayer ruthenate Sr3Ru207. It has also been 
shown that the nematic and the accompanying metamagnetic transitions are driven by an effec- 
tive momentum-dependent quadrupole-type interaction. Here, we study the microscopic origin of 
such an effective interaction. To elucidate the mechanism behind the spontaneous Fermi surface 
distortion associated with the nematic, we identify a simple tight binding model based on t2 g or- 
bit als, spin-orbit coupling and the rotation of RuC>6 octahedra as starting point, consistent with 
the Fermi surface obtained from recent angle-resolved photoemission data. Within an extended 
Hubbard model the nematic state, characterized by an anisotropy between the bands near (±tt, 0) 
and (0, ±7r), then strongly competes with ferromagnetic order but pre-empts it via a finite nearest 
neighbor interaction. We discuss experimental means to confirm our proposal. 
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Introduction — In correlated electron systems, elec- 
trons can organize themselves in states that are analo- 
gous to classical liquid crystal phases. [l| The search for 
such phases in solid-state systems, in particular for the 
quantum version of an anisotropic liquid crystal, dubbed 
electronic nematic phase, has been of great interest. Such 
a phase spontaneously breaks the point-group symmetry 
of the underlying lattice thus characteristically modify- 
ing, e.g., transport properties. Recently, a remarkable 
anisotropic longitudinal resistivity has been reported in 
the bilayer ruthenate Sr3Ru207, where the putative 
nematic phase is bounded by two consecutive metamag- 
netic transitions. [1| From a theoretical perspective, im- 
portant progress has been made as well in understanding 
the phenomenon. It has been shown that an effective 
momentum-dependent interaction model [JQ success- 
fully describes the metamagnetic transitions occurring at 
the nematic phase boundaries, Q and the behavior of the 
resistivity, when the magnetic field is applied along the c 
axis [H and tilted away from it. Q 

However, the link between the effective interaction and 
the microscopic origin is missing. While attempts to illu- 
minate the mechanism of nematic phase formation have 
been made very recently, [HI EH the focus exclusively 
lay on the quasi-one-dimensional (ID) Ru d yz and d xz 
orbitals, identifying the nematic phase as orbital order- 
ing driven by an inter-orbital Hubbard interaction. Yet, 
there is no a priori reason to exclude the two-dimensional 
(2D) d xy orbital from the picture. On the contrary, 
angle-resolved photoemission spectroscopy (ARPES) in 
the isotropic phase has revealed a van Hove singular- 
ity (vHS) near the Fermi level predominantly originating 
from the d xy orbital. [l2[ This is crucial as the isotropic 
Fermi-surface (FS) structure constrains the theoretical 
starting point for the study of the nematic phase within a 
weak-coupling theory and the nematic is associated with 
a distortion of the FS topology. [H, Q 



In this paper, we describe a microscopic route to ne- 
maticity in Sr3Ru2C>7. We first model the underlying 
band structure including all three t2g orbitals by incor- 
porating unit-cell doubling due to the rotation of the 
RuC>6 octahedra and spin-orbit (SO) interaction into a 
single layer approach. We then show that the nematic 
phase is associated with a density imbalance in the 72 
band near (±7r, 0) and (0, ±7r), where the dominant or- 
bital in the Bloch function is the d xy orbital. We find 
that nematic ordering strongly competes with ferromag- 
netic order within a multi-orbital Hubbard model, but 
pre-empts it due to a finite nearest-neighbor interaction. 
Ultimately, we discuss experimental probes to confirm 
our proposal. 

Fermi-surface topology of Sr^Ru^O-j — To obtain the 
electronic band structure we start from a single layer 
tight-binding model with three *2 g orbitals and SO cou- 
pling, as in the case of the single layer Ruthenate 
Sr 2 Ru04. 13, However, despite the similarity, two 
additional features need to be taken into account in the 
case of Sr sRi ^O?: (1) a slight rotation of the RuO octa- 
hedra [l5lll6| entailing unit-cell doubling and (2) bilayer 
coupling. Here we proceed by including (1) and discuss 
the role of bilayer coupling later. The tight-binding band 
structure H then takes the form 
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where ^ = K-sK^qM+Qs^+Q-s) in- 

sists of fermionic operators creating an electron with 
spin projection s =f, l in one of the three Ru *2g de- 
rived orbitals a = yz, xz, xy. The orbital dispersions 
are given by e£ z = —2tiCos(k x ) — 2t 2 cos(/c, / ), e^ z = 
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FIG. 1: (Color online) (a) Fermi surface of the single layer 
model. The crosses (dots) indicate the saddle points of the 72 
(71) band. The band structure parameters are ti = 0.5, t2 = 
0.05, t 3 = 0.5, U = 0.1, is = -0.03, t 6 = 0.05,^ = 0.575, g = 
0.1, A = 0.1375. (b) Fermi surface of the corresponding bi- 
layer system with finite interlayer hopping between quasi-lD 
orbitals. (c) Total DOS near the Fermi level corresponding to 
(a). See main text for details. 
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4t4COs(k x )cos(k y ) — 2t 5 [cos(2k x ) + cos(2fc w )] , 
while e]^ = —4t & s,\n{k x )sm.(k y ) describes the hopping be- 
tween the two quasi-lD orbitals and 2A J^i LA is the SO 
interaction [l7j (all energies in the following are in units 
of 2fi). For simplicity we have introduced unit-cell dou- 
bling via an effective lattice potential g with modulation 
vector Q = n). The resulting Fermi surface and den- 
sity of states (DOS) near the Fermi energy are shown in 
Figs. Ula) andQJc), respectively, where all bands are dou- 
bly degenerate due to the lack of time-reversal symmetry 
breaking. 

As shown in Fig. [T] (a), the Fermi surface is composed 
of four sheets. The hole-like a and electron-like (3 sheets 
are derived from quasi-lD orbitals, while the electron-like 
71 and hole-like 72 sheets arise from quasi-lD and 2D or- 
bital hybridization. The total DOS in Fig. [T] (c) reveals 
that the band structure provides two singularities near 
E-p, which originate from the 71 and 72 bands. Since any 
instability in a weak-coupling theory is boosted by a vHS 
near the Fermi level, these bands are most susceptible to 
a nematic transition. Consequently all other bands (a 
and 0) can be rendered less important considering the 
nematic mechanism in Sr3Ru2 07. For completeness we 
also show the effect of bilayer hopping in Fig. [T] (b). This 
FS, consistent with ARPES, exhibits a more complicated 
structure, where the CX2 sheet appears due to bilayer cou- 
pling (but not the S sheet, which most likely derives from 
e g orbitals). However, since the orbital and topographi- 
cal nature of the most singular 7 bands near (±7r, 0) and 
(0, ±7r) remains unaffected by a finite bilayer coupling, 
we proceed with the simpler single layer approach. 



Indeed, not only are the 7 bands responsible for the 
vHS but the corresponding saddle points all lie within 
shaded pocket regions I and II in Fig. Q] (a) (marked by 
crosses and dots). In region I (II), the 71 and 72 bands 
are mostly composed of only d yz (d xz ) and d xy orbitals. 
For instance, the flat parts of the 72 pockets and of the 
tips of the 71 pockets derive from the quasi-lD orbital, 
while the curved portions and the 72 saddle points have 
a strong d xy character. Note that the 71 and 72 pockets 
are formed by hybridization via SO coupling since even- 
parity perturbations such as g do not lift the degeneracy 
between quasi-lD and 2D orbitals due to different parity. 

The simple orbital composition of the 7 bands in re- 
gions I and II enables us to derive effective interactions 
within the 7 bands involving only the highlighted pocket 
regions. The singularities in the DOS suggest that ei- 
ther one or both 7 bands are involved in the formation 
of a nematic phase and that ordering can be driven by 
weak interactions due to a large DOS near the Fermi 
level. One should keep in mind, however, that there is 
more than one possible instability which takes advantage 
of the vHS. Below, we study the interactions in the 72 
band in detail and show how different orders compete. 

Microscopic route to nematicity - The bilayer 
Ruthenate has a metallic ground state. Therefore, it 
is reasonable to assume that long-range interactions are 
well screened, leaving only moderately weak on-site and 
nearest neighbor interactions. The microscopic Hamilto- 
nian then reads as 

H = H + Hi nt , 

where the extended multi-orbital Hubbard interactions 
are given by 
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with n" s = i^Ts^s t ne density operator of the a orbital on 
Ru site i. Here, U, U, and V a represent repulsive on-site 
intra-orbital, on-site inter-orbital, and nearest-neighbor 
intra-orbital interactions, respectively. To derive effec- 
tive interactions, we write H m t in the basis of the Bloch 
bands, which can be accomplished straightforwardly as 
the 7 bands in regions I and II are predominantly com- 
posed of two orbitals as discussed above. For region I, 
one thus obtains 

'akcA _ ( cos6> k -si: 
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where a\^ a and c^ a represent the annihilation of a quasi- 
particles with momentum k and pscudospin a = ± 
in the 71 and 72 band, respectively The mixing an- 
gles are given by {sin6*k, cos6>k} = 0.5 1 ± 
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FIG. 2: (Color online) Effect of an external magnetic field on 
the band structure for field orientations (a) parallel to the z 
axis and (b) parallel to the x axis with B = 0.05. 



by SO interaction. Likewise, the 7 bands can be approx- 
imated in region II by substituting the d yz for the d xz 
orbital (or, equivalently, by exchanging k x -o- k y ). 

Focusing on the 72 band, the effective interaction then 
has the form 

U\ n kl +n k2 -+ [U-2 n klCr n k2a +U 3 nk lC r«k 2 -a] 
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where n ka = c^ a c ka and the interaction strengths are 

U\ = C/(sin 2 kl sin 2 k2 + cos 2 6' kl cos 2 6' k2 ) 

+ 2[7sm 2 6> kl cos 2 6> k2 + 4V"sin 2 6» kl sin 2 k2 , 

U 2 = U(sin 2 9 kl cos 2 9 k2 + cos 2 9 kl sm 2 9 k2 ) 

+ Ucos 2 9 kl cos 2 9 k2 + 81/sin 2 6> kl sin 2 6» k2 , 

U- 3 = Usm 2 9 kl sm 2 9 k2 +U(2sm 2 9 kl cos 2 9 k2 

+ cos 2 9 kl cos 2 9 k2 ) +Wsm 2 9 kl sm 2 9 k . 2 . (4) 

Here we have included only the largest nearest-neighbor 
contribution V xy = V. 

Obviously, different instabilities compete. The dom- 
inant one is determined by the mixing parameters and 
the bare interaction strengths. To investigate the qual- 
itative features of the model we take advantage of the 
small size of regions I and II and approximate the form 
factors by their averages s 2 = (sin 2 6> k ) = 0.656 and 
c 2 = (cos 2 ^) = 0.344, which signifies the dominance 
of the d xy orbital. [l8| Hi nt then can be decoupled nat- 
urally into a nematic (n), a nematic spin-nematic (nsn), 
and a magnetic (M) channel, [l9[ 



A M = ^a^WJ), 



(5) 



with nj 11 = N 1 ^ k e j^j T n k(J and N the number of k 
points within each region. All in all, one thus arrives at 



the following mean-field Hamiltonian for the 72 band 
H MF =^Y, {E ka -V n A nPk -V nsn A nsn ap k (6) 
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where E ka denotes the 72 quasiparticle bands and p k is 
set to +(— ) for k G I (II). The effective interaction 
parameters are given by 

V n = (-0.5U + U)c 4 + Us 2 c 2 +AVs 4 , 
Vnsn = O.Wc 4 + Us 2 c 2 +AVs 4 , 
V M = O.Wc 4 + Us 2 c 2 + Us 4 . (7) 

Note that V affects only the nematic channels. While the 
above analysis is carried out for the 72 band, it is also 
valid for the tips of the 71 band inside regions I and II 
with similar effective interaction strengths. 

Effect of magnetic field and phase diagram — The 
anomalous behavior of the bilayer Ruthenate is exposed 
; in the presence of an external magnetic field, which 
can be included straightforwardly via Zccman coupling 
Hb = — B ^ i (L^+2Si). A field along the z axis, however, 
only renormalizes the spin-dependent chemical potential 
due to a quenching of orbital moments. The FS there- 
fore changes according to the dominant underlying spin 
character at each k point [see Fig. [5] (a)]. Tilting the 
magnetic field toward the x axis introduces anisotropy 
into the system via SO coupling [Fig [2] (b)]. It also in- 
creases the mixing of the r;2g orbitals in the pocket re- 
gions, where the 72 FS sheet undergoes a pronounced 
topological change and the tips of the 71 sheet are dis- 
tinctly distorted. A tilted field therefore implicates a 
low (high) longitudinal conductivity parallel (perpendic- 
ular) to the in-plane field component (not shown), in 
agreement with the experimental findings in the nematic 
phase. 

In Fig. [3] we show the mean field phase diagram as 
a function of magnetic field along the z axis and V for 
U = 0.17, U = 0.76J7, and the band-structure parame- 
ters specified in Fig. [TJ Since the rotation angles and the 
size of the RuO« octahedra vary throughout the ruthen- 



ate family, [l5j, [2JJ-|22j the nearest neighbor interaction V 



is the most eminent interaction parameter in the present 
model. As expected from Eq. iff]). V promotes ncmatic- 
ity with increasing magnitude. On the other hand, as 
V decreases, the nematic window shrinks and the phase 
boundaries, which are of first order and coincide with 
metamagnetic jumps, merge into a single line of meta- 
magnetic transitions. Driving the magnetic field strength 
up at small V/U, e.g., along the blue dashed line, tunes 
one of the 72 vHS close enough to the Fermi level to 
induce a nematic transition. This transition may also 
involve additional bands (such as 71) due to effective 
interband interactions. Note that the effect of bilayer 
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FIG. 3: (Color online) Mean field phase diagram. The order 
parameters exhibit dramatic jumps, and the magnetization 
M is shown as a red dashed line in the Am panel. Panels (a) 
- (c) display the FS evolution of the nematic 72 band near 
(±7r, 0) / (0, ±tt) along the blue dashed line. See main text for 
details. 



hopping in the presence of an in-plane field has been 
studied previously, Q but is expected to have only a 
minor influence on nematic ordering. Note also that 
in the present context nematic order characterizes the 
difference in the electronic or quasiparticle density near 
(±7t, 0) and (0, ±7r). The phenomenological order param- 
eter ~^2 k [cos(k x ) — cos(fc y )]rik used in previous studies 
|i| therefore captures the nematic phase remarkably well. 

Discussion and summary — We have described a mi- 
croscopic mechanism for nematicity in the bilayer ruthen- 
ates. Our proposal is based on a simple tight-binding 
approach including all three t 2 g orbitals, SO interaction, 
and unit-cell doubling. Leaving interactions aside, this 
approach leads to a realistic band structure with a FS 
similar to recent ARPES data. Including interactions, 
the most sensitive band is predominantly composed of 
folded d xy and unfolded quasi-lD orbitals near (±tt, 0) 
and (0, ±7r), where it undergoes pronounced anisotropic 
changes in the nematic phase. We find that the on-site 
intra-orbital interaction strongly favors ferromagnetic or- 
der over any nematic order. On-site inter-orbital interac- 
tion on the other hand slightly favors nematic order over 
both ferromagnetic and nematic spin-nematic channels 
leading to a rather non-trivial and strong competition 
between all channels. However, the nearest neighbor in- 
teraction drastically enhances charge and spin-nematic 



susceptibilities, sufficient to induce a nematic transition 
in at least one of the FS components by applying a mod- 
erate magnetic field, and thus preempting ferromagnetic 
order. 

Our results imply that nematic order is sensitive to 
the location of vHS and the balance between competing 
instabilities. The combination of SO coupling and rota- 
tion of RuOg octahedra plays a crucial role in generating 
flat bands and the vHS near the corners of the reduced 
Brillouin zone (±7T, 0) and (0, ±7r), where the band struc- 
ture is most susceptible to change. Among the members 
of the Ruthcnatc family, Ca i sS rn. 2 RuQ4 shows a simi- 
lar metamagnetic transition [24j and neutron-scattering 
pattern j2f| as non-ultra-clean S^R^Oy. 

H E3 The 

RuOe octahedra are also rotated, which distinguishes this 
compound from pure Sr2Ru04, in addition to the intro- 
duction of disorder. Therefore, it is tempting to argue 
that there also exists a hidden (due to disorder) nematic 
phase as in the case of non-ultraclean Si'3Ru207. 

Furthermore, a recent experiment probing the entropy 
landscape |23| has indicated critical behavior outside the 
nematic window. It is possible that criticality arises from 
an uncovered critical point (or weak first-order transi- 
tion) outside, yet nearby, the nematic phase space. An 
earlier work [6[ has shown that the nematic window 
shrinks as interactions decrease, and in fact becomes ex- 
ponentially small where critical fluctuations become im- 
portant. However, origin and nature of critical fluctu- 
ations outside the nematic are beyond the scope of the 
current mean-field study. 

It is worthwhile to comment on the difference between 
our proposal and the one in Refs. [l(| and [llj], where 
nematic order is associated with orbital ordering between 
the two quasi-lD orbitals, and the anisotropy appears in 
the cti and «2 FS sheets near the T point. A possible 
way to distinguish both proposals is scanning tunneling 
microscopy, [28|, [2^] which in principle should probe all 
bands. Quantum oscillations would be another way [30I ] 
since the 72 frequency should drastically change between 
the two isotropic phases indicated in Figs. EJa) andEJc), 
and, although difficult to observe, should split into two 
when passing through the nematic phase. 
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